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FOREWORD 


This report covers part of the work performed during the second year 
phase of a two year research project, at the School of Mechanical Engineer- 
ing at the Georgia Institute of Technology in Atlanta, Georgia, for the 
NASA Manned Spacecraft Center, Houston, Texas. The contract designation 
is NAS 9-10415 and the project title is "Study of Design Parameters of 
Space Base and Shuttle Heat Rejection Systems." The project was monitored 
by Dr. W.E. Simon of the Power Generation Branch of NASA MSC, Houston, 
Texas, and was performed by Dr. W.Z. Black and Dr. W. Wulff as Co-Invest- 
igators. The project resulted in one Annual Report [2]* and two Final 
Reports, the first one of which covers a detailed and rigorous space 
radiator simulation analysis [1] and includes a computer program users 
manual. The second Final Report is this report. 

Although the study of system parameters and the system optimization 
were originally conceived to be performed on the basis of the rigorous 
radiator system simulation [1], the growing complexity of this simulation 
soon gave rise to the need for a simplified analysis. The simplified 
system simulation was later expanded into systematic optimization pro- 
cedures. Both the simplified simulation and the optimization procedures 
serve to supplement and support the originally intended system parameter 
study [1], 

The work presented here was supported by the contributions of computer 
coding by Mr. Richard J. Huntley and Mr. Wallace W. Carr, both Graduate 
Research Assistants and M.S. Candidates. 


*Numbers in brackets refer to the Bibliography at the end of this report. 
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SUMMARY 

A simplified radiator system analysis was performed to predict steady-state 
radiator system performance. The system performance was found to be describable 
in terms of five non-dimensional system parameters. The governing differential 
equations are integrated numerically to yield the enthalpy rejection for the 
coolant fluid. 

The simplified analysis was extended to produce firstly the derivatives of 
the coolant exit temperature with respect to the governing system parameters and 
secondly a procedure to find the optimum set of system parameters which yields 
the lowest possible coolant exit temperature for either a given projected area 
or a given total mass. The process can be inverted to yield either the minimum 
area or the minimum mass, together with the optimum geometry, for a specified 
heat rejection rate, 

'The major accomplishments of the simplified radiator system analysis are: 

(1) the reduction of the number of necessary systems parameters 
from twelve or more to six, 

(2) the graphical representation of system performance in terms 
of non-dimensional groups, suitable to aid in the design of 
radiative heat rejection systems, and 

(3) an efficient computer code suitable for preliminary 
performance prediciton. 

The accomplishments of the systematic optimization analysis are two com- 
puter codes which perform iterative optimizations processes leading to the 
maximum heat rejection for either a given projected fin-plus-tube area or a 
given total system mass. 
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I . INTRODUCTION 


A large-scale, complete and rigorous computer simulation of space 
radiator systems was developed under the same contract as the work 
presented here [1],* The computer code consists of over fifty program 
units and is capable of simulating transient as well as steady-state 
radiator system performances under prescribed time-dependent operational 
and environmental conditions. The program accommodates both gaseous and 
liquid coolant fluids with any consistently prescribed set of thermodynamic 
and transport properties. In principle, this large-scale computer program 
could serve to not only predict radiator system performance but also to 
optimize certain design paramters via enumeration of performance charac- 
teristics, associated with selected parameter combinations 

Any large computer simulation must be tested during its development; 
the greater the number of independent verification modes the greater will 
be the confidence in the program performance. In addition to the tests 
described in Reference [1], a simplified radiator system analysis was 
developed to verify the large-scale computer program performance. This 
simplified analysis served later (1) to extend the rigorous analysis to 
a wider class of radiator geometries and operating conditions [1] than 
would have been possible otherwise and (2) to develop systematic optimi- 
zation procedures in support of system parameter studies on the basis 
of the rigorous computer simulation. 

The essence of the simplified analysis lies in the recognition of 
the dominant performance characteristics and in the parameter reduction 
through normalization of the governing differential equations. The 
optimization is based on standard requirements of extrema subject to 
suitable constraints. These requirements are imposed on the governing 
differential equations prior to their numerical integration. 

Chapter II below is a discussion of the simplified analysis and is 
followed by the presentation of the optimization analysis in Chapter III. 


^Numbers in brackets refer to the Bibliography. 
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II. SIMPLIFIED RADIATOR SYSTEM ANALYSIS 
Ao Objective and Background 

The purpose of this analysis is to provide an efficient process of 
least complexity which serves to describe the steady-state performance of a 
space radiator system. The radiator system is described in detail in the 
following chapter. 

The analysis is intended to accommodate the essential features of 
the radiator system performance but to require less computational effort than 
the rigorous system simulation described in Reference [1] . Finally, the 
simplified analysis should serve as the basis for systematic optimization of 
the radiator system. 

A first simplified analysis was developed during the first year of 
the contract period [2], also in support of the rigorous computer analysis, 
and for the purpose of treating approximately some radiator system geometries 
which deviate from the basic geometry underlying the rigorous system analysis, 
such as cylindrical shell radiator structures, asymmetrically loaded coolant 
channels etc. [l] . This analysis was based on two global energy balances, 
one for the coolant fluid and one for the radiator panel, and included the 
two thermal resistances associated with the energy transfer from the fluid to 
the tube and from the radiator into space. The mean tube wall temperature 
(and with it, the fin base temperature) was obtained from a single, transcen- 
dental equation (Eq. D.ll, [2]) and the cooling capacity was found explicitly 
in terms of the mean tube wall temperature (Eq. 1.1 or 3, [2]). Agreement 
between this simplified analysis and the full-scale computer simulation was 
found to be approximately 47 c . 

The search for an optimum radiator geometry on the basis of the 
first simplified analysis lead to the conclusion that for optimum radiator 
geometries neighboring coolant channels should be expected to be close 
enough to each other so as to obstruct each other’s, and the fin panel’s, view 
of space. Moreover, the rigorous computer simulation indicated that the 
tube wall temperature changes significantly in the direction of the flow even 
though conduction parallel to the flow direction, both in the fluid, in the 
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tube wall and in the fin, is negligible. These effects are not accounted for 
in the first simplified analysis. 

The new simplified analysis presented here, for the prediction of 
steady-state radiator heat rejection performance, takes into account 

(i) temperature variations within the fluid, tube 
and fin, along the direction of the flow, 

(ii) temperature variation within the fin, normal 
to the direction of the flow, 

(iii) direct radiative interaction between fin panel 
and space, tube and space, fin and tubes, 
tube and neighboring tubes (by approxima- 
tion) , 

but does not take into account 

(i) thermal conduction in the direction of the 

coolant flow, both in the fin-tube structure and the fluid 

(ii) temperature variation of thermal properties, 

(iii) end effects at the coolant fluid manifold, 

(iv) cross-sectional changes in the fin panel, 

(v) secondary reflections of radiant heat. 

The analysis reduces to an initial value problem consisting of two 
ordinary, first-order, coupled differential equations, with one of its initial 
conditions given through a transcendental equation. The new formulation is 
suitable for optimization,. 
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B. Description of Radiator System 

The radiator system considered here consists of a given number of 
parallel coolant channels, equally spaced in one plane, and connected by 
rectangular fin panels which have constant thickness and are symmetric with 
respect to the plane through the tube areas. 

The coolant channels are taken to be thin-walled tubes of circular 
cross-section, with diameter d and length L. Substitution of the hydraulic 
diameter for the diameter permits readily the inclusion of channels with non- 
circular cross-sections, except for possible effects from radiative inter- 
action between tube, fin and environment. 


Let the distance between tube centers be designated by d + 2H and 
the fin panel thickness by t. Then the radiator can be considered to con- 
sist of the given number n of identical fin elements as shown in Figure 1. 


The coolant fluid enters the channel at z = 0 with time-invariant 
temperature T^ = T q . Both the fluid temperature T^(z) and the fin base 
temperature T^Cz) decrease along the tube provided the equivalent sink tem- 
perature T g is less than T q . The coolant fluid emerges from the radiator at 
z = L with the exit temperature T^. 

The objective of the analysis is to predict the rate of heat 
rejection m (T q - T^) per tube for a given mass flow rate m per tube and 
given fluid properties, namely density p, specific heat c , thermal conduc- 
tivity k^ ? dynamic viscosity U, and given properties of the fin, that is, 
surface emittance £, solar absorptance 9 and thermal conductivity k. 
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Coolant Exit 



Figure 1. Geometry of Fin System Element 
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C . Governing Equations 
1. Conservation of Energy 

Setting the local change of the fluid enthalpy flux first equal 
to the local heat exchange per unit tube length between fluid and tube wall, 
and then equal to the radiative heat rejection per unit tube length one ob- 
tains, for i effective sides of the radiator: 


. dT f 

" m C P dT = " dh c (T f - V 


(i) 


ndh (T, - T, ) 
c r b 


2iHean(T^ - T 4 ) 
b s 


( 2 ) 


where h^ and g~ represent, respectively, the convective film coefficient and 
the S tefan-Boltzmann constant. The symbol rj represents the overall tube- 
and-fin effectiveness and is evaluated in Section 3. The convective film 
coefficient is computed from the Nusselt number N 

Nu 


h 

c 


= K 

Nu d 


(3) 


the evaluation of which is deferred to Section 2. The equivalent sink tem- 
perature T g is computed from the known incident normal radiant heat flux q 11 
by 


T 

s 



(4) 


The first-order differential equation, Eq. 1, is subject to the initial con- 
dition at z = 0 


T-(0) = T 
r o 


( 5 ) 
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Equations 1 and 2 determine the two dependent variables T^(z) and T^Cz)- 
Integration of Eq. 1 from z = 0 to z = L gives the unknown fluid exit tem- 
perature T g . 
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2. Convective Film Coefficient 

The Nusselt number N^ u in Eq. 3 depends on the Reynolds number 

N , the Prandtl number and the d/L ratio. 

Re rr 


N 


Re 


4m 

irdy 


( 6 ) 


N 


Pr 



(7) 


For laminar flow, N Rg ^ 2300, of non-metallic fluids, Np^ ^ 0.1, 
Hausen [3] established the relation 


0.0668 N N - l 
Re Pr L 


N, 


Nu 


3.65 + 


2/3 


1 + 0.045(N Re N p -) 



0.14 


( 8 ) 


The nomogram in Figure 2, taken from the VDI Waermeatlas * facilitates the 
estimate of the Nusselt number in accordance with Eq. 8* 

For turbulent flow, N > 2300, of non-metallic fluids, N_ ^ 0.1, 

5 Re Pr 

similarity considerations lead to 



which is also presented in a nomogram in Figure 3. 

Liquid metal convective heat transfer in tubes, Np^ <0.1, may be 
represented by the result of the work by Seban and Shimazaki [4] which is 


*VDI-Verlag GMBH, Dusseldorf, West Germany 
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valid for (N_ N„ ) > 100 and L/d > 60 

Re Pr 


N._ =5.0+0.025 (N N ) 0,8 

Nu ' Re Pr 


( 10 ) 


Other Nusselt number relationships may be more suitable in a 
particular case. Their selection does not affect the ensuing analysis be- 
cause the Nusselt number calculation is part of the input data preparation, 
prior to the numerical evaluation of the solution. 
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Figure 2. Nussett Number Nomogram for Laminar Flow of Non-Metallic 
Coolant Fluids, Eq. 8. (VDI Waermeatlas) 
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Figure 3. Nussett Number Nomogram for Turbulent 'Flow of Non-Metallic 
Coolant Fluids, Eq. 9. (VDr'.Waermeatlas:).: 
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3. Total Fin-Tube Effectiveness 

The symbol n in Eq. 2 denotes the overall radiative fin-and-tube 
effectiveness and accounts not only for the non-uniform temperature dis- 
tribution within the fin and normal to the flow direction, but also for the 
effects on the radiative heat transfer from the tube. These effects are 
due to the additional radiating tube area and the partial blockage of emis- 
sion from the fin and the tubes by the tubes of the neighboring fin elements. 
Adding the contributions from the tube to that of the fin, one obtains 

n(z) = M(z) ’n(z) + N (11) 


Here, n(z) is the radiative fin effectiveness of the unobstructed fin* Eqs. 32, 
[5,6,7] and varies with the fin base temperature T^(z). The factor M is 
the product of the shape factor of the fin with respect to the sink, F^, 
and the corrective term f on n for the temperature profile distortion re- 
sulting from the tube-fin interaction [8] : 


M = F fg (X)*f(X,N c ) 


where 


( 12 ) 


X 


4H 

d 


(13) 


and 


N 

c 


2eoH 2 

kt 


[T b (z)] 3 


( 14 ) 


which is the well-known conduction parameter. Hottel’s crossed-string 
method yields 
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F fs (X) - 


U 


A (A + 2) + arc sin --"y * ~ */2 I /A 


(15) 


and the curve fit of data in Reference [8] resulted in 


f (A ,N ) = 1 - (N /A) (0.1460 N - 0.02866) . (16) 

c c c 


Finally, the symbol N in Eq. 11 represents the radiative heat re- 
jection contributed by the tubes and is the shape factor of the tube with 
respect to the sink, multiplied by the tube-to-fin area ratio 


N(A) 


= it 1/2 + [ A + 2 - y A (X 


+ 2) - arc sin 


■it]'"] 


/A 


(17) 


This completes the description of the overall fin-and-tube 
effectiveness. In summary it should be emphasized that axial conduction, 
radial tube wall temperature variation, kinetic and potential coolant 
energies, and end effects are ignored, that radiative interaction between 
tube and fin and between tubes is only approximate and that the represen- 
tation of Eq. 16 is limited to T < 0c8 T, (but could be extended, in 

S D 

principle) . 
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D. Scaling Parameters 

Effective analysis, numerical integration, graphical representa- 
tion of numerical results and, ultimately, the system optimization makes 
the reduction of parameters mandatory. Introduce the non-dimensional 


axial distance 


temperature 


reference conduction 
parameter 


ratio of convective to 
radiative resistance 


convection number 


C " L 

(18) 

0 = T/T 

o 

(19) 

3 2 

2eaT 0 H 
kt 

(20) 

itk 

' * “ k £ N NU 

(21) 

■ *V N G z 

(22) 


where the Graetz number N 

Gz 


is defined by 


N = (c m) / (k L) 
b z p r 


— — N N 
4 L Re Pr 


(23) 


With these parameters and the A defined by Eq- 13 one may recast the problem 
as previously established by Eqs„ 1,2 and 5 as given by Eqs„ 24 through 26: 


d0 

3T ' - u(e f ■ V 


(24) 


6 f - 9 b * V "< 6 b - 6 J> 


(25) 



15 


subject to the initial condition at C = 0 


6 f C0) = 1 . (.26) 

Integration of Eq. 24 subject to Eqs. 25 and 26 yields the coolant fluid exit 
temperature 9 g sought 


e 

e 


= e f (i) 


(27) 


and, consequently, this coolant fluid exit temperature is a function of five 
parameters 


0 

e 


- 0 e (U,V,N c ,A;0 s ) 


(28) 


The first four parameters correspond to the geometrical dimensions L, H, t 
and d of the radiator system and 0^ represents the environment. 
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E. Solution 

Equation 25 is transcendental because of n(0^) and would require 
iterative solutions at every step of the numerical integration of Eq. 24. 
It is more economical to derive a second differential equation for 6^ by 
differentiating Eq. 25 


de u 

b 

dC 


- U 




+ 



(29) 


and to obtain the initial condition for Eq e 29 through a single iterative 
process (Newton Raphson method) from 


Vn(e - e ) + e = l at 5 = o 


(30) 


The derivative db/dO^ is obtained from Eqs c 11 and 16 


dN 

dn _ dn c 

de, ~ dN d0 u 
b c b 


3 N 0 K 2 (Mn) 
cb 3N 

c 


= 3N 0, 
c b 


nF 


fs 


0 . 


02866 - 0.2920 


N 0 ' 
c b 


/A 


(31) 


+ M 


dn ( 

dN c I 


which can be evaluated once the fin effectiveness o(N c ), that is the ratio 
of the actual power loss from an unobstructed (.tubless) . fin panel to the 
power loss from an ideal, unobstructed fin of infinite thermal conductivity, 
is represented by a power polynomial: 



17 


6 

n(N c ) - a^N/ 

i = 0 


* (32a) 


where a 0 = 1.000 000 

a 1 = - 1.163 143 

a 2 = 1.478 836 

a 3 = - 1.267 550 

a 4 = 0.632 522 3 

a 5 = - 0.162 706 7 

a 6 = 0.016 542 23 

Equation 32 is valid for 0 ^ N c £ 2.5. For greater values, N £ > 2.5 

- 0.229 771 8 

n(N c ) = 0.686 609 5 e (32b) 

Equations 24 and 29 have been integrated by a fourth-order Runge- 
Kutta procedure [2], subject to initial conditions given by Eqs. 26 and 30. 
The result of the integration yields primarily 0^ in the form of Eq. 28. 
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F. Results 


A typical computer print-out of the results is shown in Figure 4. 
First are presented the five parameters { 0 ; U,V,N , A} of Eq. 28, as read- 
in. Following that is a table in which are listed, as functions of axial 
distance measured from the coolant fluid inlet, the fluid temperature 0^ = 
THETAF, the fin base temperature 0^ = THETAB, next the local conduction 
parameter = NC, the local total fin-plus- tube effectiveness n = ETABAR, ■ 
the local fin effectiveness n = ETA, the nondimensional axial position £ = 
ZETA and, finally, the number of integration steps required to reach the 
particular axial position. The number of integration steps depends on the 
chosen error limits and the rates of variable changes. 

The time rate of heat rejection per fin element (tube) is to be 
computed from 


me (T - T ) = m c T (1-0) 
p o e p o e 


(32) 


Parameter studies are presented in Figures 5 through 9. The 
non-dimensional fluid exit temperature is plotted versus non-dimensional fin 
panel area with the temperature 0 g (Fig, 5), the non-dimensional tube 
spacing X (Figs. 6 and 7) and the reference conduction parameter (Figs. 

8 and 9) as parameters- 






a 



0 ajn^Baaduiax 3TX3 }ubxoo3 x buo T suauii g-uo^ 


Non-Dimensional Fin Area A 

Figure 5. Non-Dimensional Coolant Exit Temperature 0 vs, Non-Dimensional Fin Panel Area 
with Sink Temperature 0 As Parameter. 
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Figure 9. Non-Dimensional Coolant Exit Temperature VS. Non-Dimensional Fin Area Fi 
Varying Conduction Parameter and Short Fins 
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III. OPTIMIZATION OF THE RADIATOR SYSTEM 
A. Purpose 

In the design of the radiator system all parameters describing the 
system objectives and the system environment are prescribed beforehand while, 
at some occasions and within certain limits, the geometry of the system may 
be selected in the process of the design. This possible freedom of choice 
leads to the optimization of the radiator system geometry. 

The primary purpose of the radiator system is to lower the coolant 
fluid temperature. For a given set of constraints the radiator with the 
lowest fluid exit temperature is considered here to be an optimum radiator 
system. 

The set of constraints to be imposed on the optimum radiator system 
must be selected for a particular set of circumstances. In this analysis 
two sets out of many possible constraints were selected, that is, weight 
(or volume) and projected area of the radiator are considered to be prescribed. 

The process of optimization proposed here, namely to find the 
minimum coolant fluid exit temperature for either a fixed radiator area or a 
fixed m??s may be inverted to yield the minimum radiator area or weight for 
a prescribed cooling rate. 

The optimization process can be outlined in the following steps: 

(i) select, for given coolant exit temperature 

e e = 0 e {u , v ,N c ,X; 6 s > (28) 

that optimum set {U,V,N c ,A} of geometric parameters which 
constitutes, under selected constraints, the system geometry 
for the lowest possible coolant exit temperature, 

(ii) transform the non-dimensional groups U, V and N into suitable. 

c 

non-dimensional radiator length L, fin height H* and panel 
thickness t,* 
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(iii) calculate the physical dimensions L, H, t and d of the 
optimum radiator system. 

If all fin elements operate under optimum conditions then the 
radiator as a whole operates under optimum conditions. 
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B. Problem Formulation 

1. General 

Of the five parameters U, V, N , A and 0 which define 0 only 

c s e 

the first four can be considered independent while 0 g is determined by the 
incident radiant flux. For simplicity let us introduce the vector X with 
components x^ = U, X£ = V, x^ = and x^ = A. The minimum coolant fluid 
exit temperature subject to r < 3 constraints may be obtained from 

the necessary conditions for a relative extremum 

00 

= 0 i = 1, .... (4 - r) (33) 

i 

and the r constraints 

cPj(X) = 0 j * i, ...» r . (34) 

This problem formulation implies that the r constraint equations can be 
solved explicitly for r components of X which can be substituted into 
Eqs. 33 to yield (4 - r) equations for the remaining (4 - r) components. 

This was in effect achieved by replacing the Nusselt number relationships 
given in Eqs. 8 and 9 by a simpler relation of the form 

\u - c V “pr" (d/L)P (35) 

with appropriate constants c = 0.023, m = 0.8, n = 1/3, p = 0, for >_ 2,300; 
and c = 1.86,m = n = p = 1/3 for N R ^ < 2,300. The method of Lagrangian 
multipliers, however, permits the retention of the original constraints at 
the expense of additional, computational efforts in the case of fixed-mass 
constraints, since Eqs. 34 need not to be solved explicitely for X. 

Two sets of constraints have been considered. Both sets have in 
common four sets of bounds 
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U . <_ x < U 

mm 1 max 


V . <. x < V 

mm 2 max 


(N ) . < x < (N ) 

c min 3 c max 


* . < x < X 

mm 4 max 


(36) 


which represent practical limits of the design parameters x^, i = 1, •..,4. 
Should the search procedure lead to the intersection of any of the above 
limits, then their appropriate equalities would enter the set of constraints. 
The other two constraints, namely of fixed volume or of fixed area, are 
discussed below. 
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2, Area Constraint 

Introduce the characteristic length of a single fin element 

(37) 
o 

which may be interpreted as the length of a square fin element 
whose effectiveness is unity and which reduces the fluid temperature to zero. 
Dividing the physical dimensions L, H, t and d by the reference length L q one 
obtains the non-dimensional geometric quantities L, H, t and d, respectively. 
The normalized projected fin element area is (for t <<d) 

A* = A/L 2 = [ 2HL + dL]/L 2 . (38) 

o o 

Consequently, the constraint of "fixed area," the first of Eq. 34, becomes 

cp 1 = 2H*L*[ 1 + 2/ A] - A* = 0 (39) 

or 

cp 1 (x 1 ,x 2 ,x^) = XjX^l + 2/x^) - A* = 0 (40) 

This constraint equation can be solved explicitly for any one of its arguments, 
regardless of applicable Nusselt number relationship (see comments following 
Eq. 34). 

The constant area constraint may be supplemented by an additional 
constraint reflecting the preference for a particular fin panel thickness, 
arising from manufacturing considerations. A fixed panel thickness t, or 
equivalently t" would require 



cp 2 (X) = t* - ♦ • N 2 u (X) = 0 


(41) 
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where <j> 


becomes 


is known beforehand 


= (1/ 2i) (Ttk^/k) 2 k 2 / (im C e a T 3 ) 

l t p o 


(42) 


Equation 41 depends on the Nusselt number relationship used. It 
a necessary constraint when the optimization tends toward panel thick- 


nesses which approach the tube diameter. 
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3. Volume Constraint 

The volume constraint is stipulated by the desire to design a 
radiator system of minimum mass. The coolant mass may (liquids) or may not 
(gases) contribute significantly to the total mass. 

The mass of one fin element is, under the assumption of a tube 
wall thickness, equal to the fin panel thickness 


* 

M = 


M 


+l M o 


* * 
L t 


it it it 

[2H + ir(d -t )] = 




(43) 


provided t << d . Here M o =pL Q is the reference mass and M the physical 
mass of the structure. With 


<f> 2 = 2i(p f /p) (k/k f ) 


the coolant fluid mass may be expressed as 


it 



M, 


♦l M o 




U 


(44) 


(45) 


Consequently, the "volume constraint" becomes 

*!<?> = <ot “ *L X 2 N Nu ( ? } j X 2 U + 27r/x 4 ]N Nu /x 3 + *2 /x 4 j = ° (46) 

it it it 

where M = M + M^. Equation 46 reveals that while ^ > 1 the fluid mass 
contributes to the total mass only when x^ becomes small. For the expected 
cases of x^ * 10 Eq. 46 reduces to a volume constraint as the influence of 
the density ratio p/p^ vanishes. 
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where 6 is a small vector, with components {Ax^ Ax^}, an ^ with its 

endpoint on the hypersurface defined by Eqs. 34. 

There exist analytical expressions for the sufficiency conditions 

2 

in terms of second-order derivatives 3 6 e /3x^3Xj which are applicable when 
either all constraints can be solved explicitly as implied by Eqs. 33 and 
34 or the number of constraints is at least two. Since, however, 0 is 
obtained through numerical integration of Eqs. 24 and 29 and the sufficiency 
test is to be performed only once the potential optimum is found, it appears 
economical to evaluate Eq. 47 directly. Further developments concerning the 
sufficiency criteria are necessary at this time. 


3 
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C. Solution 


1, The Optimum 

The remaining task is to solve the system of Eqs. 33 as they are 
obtained after substitution of Eqs. 34. The solution is obtained through an 
iterative process based on the Newton-Raphson procedure. Starting with an 
estimated set of parameters the iteration is carried on according to 

\ , k = 1,2, ... (48) 


where represents the current, X^ + ^ the future parameter set and the 
increments AX^ are the solution to the system of linear algebraic equations 


\ - <*> k <4 ?>k • 


(49) 


The current components (y.) of Y, are the current values of the derivatives 

1 k 


30 
e 

3x. 


, i = 1, . . . , (4 - r) 


(50) 


The current elements {A..} of the square matrix (A), are the 

ij k * k 

second-order derivatives 


A. . 


2 

3 9 


e 


3x^3x^ 


i ,3 = 1, 


(4 - r) . 


(51) 


It is obvious that the matrix A^ is symmetric and one needs to compute only 
(5 - r) (4 - r)/2 independent, second-order derivatives. 

The derivatives in Eq. 50 and 51 are obtained from Eqs. 24 and 25 
by first differentiating with respect to x^, i = 1, ..., (4 - r) and then 
interchanging the order of differentiation. This leads first to 2(4 - r) 
first-order, ordinary, non-linear differential equations for 
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d_ 

d? 




, i = 1, ...» (4 - r) 


, i = (5 - r), (8 - 2r) 


(52) 


and then to (5 - r) (4 - r) first-order, ordinary non-linear differential 
equations for 


dA. . 

d 

/ A f \ 

d? 

d? 

l s V x j ) 

fn 

= <L 1 

{ ^ ) 

d4 

dC 

\ s V x j / 


I (53) 


The second sets of equations in Eqs . 52 and 53 are necessary 
because of the dependence of 0 g on 0^. All initial conditions, at £ = 0, 
for Eqs. 52 and 53 can be derived from Eqs. 26 and 30 


y^o) = (y D ) 


V 0) = (A o } ij 


, B..(0) = (B o ) 


(54) 


ij 


In summary, one needs to integrate, together with Eqs. 24 and 29, 
the Eqs. 52 and 53, subject to the initial conditions given by Eqs. 26, 30 
and 54. The integration is carried out from £ = 0 to £ = 1 where it yields 
not only 0^ and 0^(1) but also all derivatives in Eqs. 49. This system is 
solved for the components of AX, then a new set of parameters X is computed 
from Eqs. 48 and the iterative cycle repeated. The repetitions are continued 
until either the first (4 - r) y^ f s are small or one of the inequalities in 
Eqs. 36 is violated. Should that happen then an additional constraint is intro- 
duced, r is incremented by one, and the cycle is continued until either all 

limits in Eqs. 36 are reached or the remaining y. f s are sufficiently small. 

1 * * 

The result may or may not be an optimum for the initially chosen A or M . 

Finally, the potential optimum is tested in accordance with Eq. 47. 
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2. Parameter Transformation and Physical Dimension of Optimum Radiator System . 

Once the optimum set -£x^ , . . . ,x^ is found it is necessary to compute 
from this the physical dimensions L, H, t and d which define the geometry of 
the radiator system. To accomplish this task the set jx^ , . . . ,x^ first 

transformed into the previously introduced set of normalized dimensions 
•|l*, H*, t*, d*J. This transformation is, in general, not possible in explicit 
form. However, when the simplified Nusselt number relation of the form given 
by Eq. 35 is substituted for the more general expressions in Eqs . 8 and 9 then 
L*, H*, t* and d* can be expressed explicitely in terms of the parameters 

>••• > x 4 J* 

Define first the starred quantities. 



c (N* e )m(N pr )n 


(.57) 



(58) 


which are essentially the unstarred quantities evaluated with the reference 
length instead of any other dimension and which are all known from system 
objectives and environmental conditions. 

Next one needs to distinguish be tween laminar and turbulent coolant 
flow as the constants c, m, n and p are different for the two regimes. 

F or laminar flow, c = 1.86, m = n = p = 1/3 



( 59 ) 
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H 


* 



(60) 


t 


* 


tt 2 cN*/(N * z ) 2 



X 1 X 3 


and 


(61) 


d 


* 



* 

Gz 



(62) 


For turbulent flow, c = 0.023, m « 4/5, n = 1/3, p = 0. 

With 


*3 ’ ’ N Nu /B Gz 

(63) 

obtains 


M 4/9 


L* - 

(64) 

H* = 1/2 [2 5 4 > 3 14 ] 1/9 (x 2 5 x 4 4 ) 1 / 9 

(65) 

(X 5 x 4 ) 2 / 9 

t* = N*/4[2 5 (() 0 14 ] 2 / 9 — - — 

c T 3 x^ 

(66) 

* r x ?i 5/9 



(67) 

Finally, multiplication of the starred quantities L*, 

H*, t* and d* 


by the reference length gives the dimension L, H, t and d of the fin 
radiator system with optimum performance. 
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D. Results 

The ultimate presentation of the results from the radiator system 
optimization is the plot of the optimum non-dimensional geometric para- 
meters L*, H*, t*, and d* versus the non-dimensional enthalpy rejection 
(1 - 0 e ). These graphs would include three or more parameters, namely 
the sink temperature and the parameters occurring in the constraint 
equations. However, the optimum geometry is not expected to depend 
strongly on all parameters because of physical considerations, and con- 
venient graphical presentation should be possible. 

As indicated in the Introduction, the systematic optimization based 
on analytical extrema search techniques was developed in addition to the 
original program objectives but could not be completed within the contract 
period. The graphical presentation of the optimum geometry in terms of 
intended enthalpy rejection is therefore not included in this report. 

However, computer codes were written on the basis of the solution 
discussed in Section C. These codes produce the optimum radiator para- 
meters L*, H*, t* and d* based on the Nusselt number relationship in Eq. 
35, for any particular set of input parameters computed from mission re- 
quirements and environmental conditions. These results are described in 
detail in Chapter IV, Sections C. 4 and D. 4. 
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IV- COMPUTER CODES 
A. Introduction 

Three separate F0RTRAN codes were developed, one for the simplified 
radiator system simulation, one for the minimum area optimization and one for 
the minimum mass system optimization. 

The three codes have several subprograms in common and could be 
united into a single code to avoid duplication. The simulation program, 
however, is more efficient as a single program and will remain a tool by 
itself. The optimization programs, on the other hand, are means to develop 
a tool, namely suitable working charts from which to read optimum system 
parameters. Once these diagrams are obtained, the codes are no longer needed. 

The three F0RTRAN codes are discussed separately in the following 
three sections but reference is made in the description of subprograms 
which have previously been discussed. 

B. Radiator System Simulation 

1. Objective 

The purpose of the Simplified Radiator System Simulation (SRSS) 
is to evaluate the analysis developed in Chapter II for determining the coolant 
exit temperature 

0 e (U,V,N c ,X;O g ) 

in terms of the five governing groups U, V, N , A and 0 defined by Eqs. 13, 

c s 

19 through 22. 

The coolant exit temperature and also the fin base temperature 
are obtained by integrating simultaneously Eqs. 24 and 29, subject to the 
initial conditions given by Eqs. 27 and 30. The integration is performed 
by the Runge Kutta integration procedure described in Reference [1], under 
variable step size mode. 
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The individual program units are discussed in the subsequent 
sections to the extent deemed necessary for the proper utilization of the 
code. Following the program description are presented the input and output 
specifications . 

2. Deck Assembly 

The Simplified Radiator System Simulation (SRSS) code consists of: 


one 


main program 

MAIN 

three 

SUBR0UTINE 

subprograms 

RKS 

SDRV 

SCNTL 

three 

FUNCTI0N 

subprograms 

ETA 

P0LY 


A block diagram is shown in Figure 7 below. 

a. The MAIN Program accepts data input and lists the accepted input 
data. Following the input data management, the initial conditions for inte- 
gration are set, an initial step size is computed and control variables 
are set to control the integration procedure. Then the Runge Kutta procedure 
RKS is called which performs the integration and indirectly the listing of 
results. . After return from the integration, control is transferred to the 
start of the program for arbitrarily many repetitions of the program execution. 
When all input data are exhausted then control is transferred to a normal 
exit. 

The input data preparation is discussed in Section 3 below. 

The initial conditions are given through Eqs. 27 and 30. The 
dependent variables 0^ and 0^ are, during integration, placed in the array 
Y; Y(l) = 0£ and Y(2) = 0^. Equation 27 leads to the statement Y(l) = 1.0 
and Eq. 30 is solved for 0^ by the Newton-Raphson iteration which follows 
the listing of accepted input data. The iteration process leads to the 
statement Y(2) = THETAB (= 0^ ) . 




Figure 10. Block Diagram For Simplified 


Radiator System Simulation Code 
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Control variables which direct the Runge Kutta Integration 

RKS are: 


zzz 

= 

0.0 

DX 

= 

1/(U(1 - 

IFVD 

= 

0 

N 

= 

2 

IBKP 

= 

1 

NTRY 

= 

1 

IERR 

= - 

0 

A(I) 

= 

5 x 10 -5 

R(I) 


5 x 10 -5 


initial value for independent variable 

initial step size 

variable step size selection 

two equations to be integrated 

cut interval as required 

normal integration mode 

normal error condition (DX ^ 0) 

absolute error per step 

relative error per step 


are: 


Control variables which direct output during the integration 


DLMT = DXWRT/500 
XWRT = 0.0 

LSTEP = 0 
LCNT = 0 
ICNT = 0 


error limit for output interval 

first value of independent variable at 
first output listing 

controls change to fixed step, LSTEP = 1 
once integration step exceeds write 
interval DXWRT 

senses cut-back (LCNT = 1) of step size 
prior to listing of results; used to 
save uncut time step 

integration step counter 


The CALL statement for the Runge-Kutta SUB-R0UTINE is explained 
in Reference [1] together with the specification of necessary DIMENSI0N and 
EXTERNAL declarations. SDRV and SCNTL must be declared EXTERNAL. The 

necessary array allocations can be read from the listing of the RKS code: 
all arrays with variable dimension (N) must be declared in MAIN for N = 2 
elements. Only the elements of Y, A and R must be specified prior to the 
CALL RKS statement. 
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b. The Runge-Kutta SUBROUTINE RKS is described in detail in Reference 
[l]. The user is not expected to alter this code in any way and needs to 
know only the control features discussed in paragraph a above and the fact that 
RKS calls two subroutines SDRV and SCNTL. 


c. The SUBR0UTINE SDRV serves to evaluate the derivatives d0^/dC = dY(l) 

and d0, /dr = dY(2) in accordance to Eqs. 24 and 29. The fin effectiveness 
b 

T) (N ) and its derivative dT/dN are evaluated externally by the FUNCTI0N 
c c 

subprograms ETA and DETA. 

SDRV is called from RKS and supplies the derivatives, via argument 
list, as functions of the independent variable X and dependent variables Y(l) 
and Y(2), to the integrating procedure RKS. SDRV is called four times per 
integration step. 

d. The SUBR0UTINE SCNTL serves to provide 


(i) output listing at prescribed intervals DXWRT, 

(ii) transfer to fixed-step integration mode once 
the automatically selected (by RKS) step size 
DX exceeds the print-out interval DXWRT, 

(iii) termination of the integration procedure and RETURN 
to the MAIN program (via RKS). The same program, 
with different WRITE commands is used in the 
Minimum Area and Minimum Mass Optimization codes. 


In the argument list are available 


(i) 

(ii) 

(iii) 

(iv) 

(v) 


the dependent variables Y(l) and Y(2) 
their derivatives DY(1) and DY(2) 
the dependent variable X 

the step size DX which may be altered for 
repeat of last step 
NTRY = 1 normal continuation of RKS 
2 return from RKS to MAIN 
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(vi) 


3 

4 

IFVD = 0 
1 


repeat with new DX 
restart 

variable integration step 
fixed integration step 


e. The FUNCTION ETA Subprogram computes the fin effectiveness T|(N^ ) 
for the unobstructed fin in accordance to Eqs. 32a and b. The power poly- 
nomial is evaluated by calling the FUNCTI0N P0LY subprogram [1]. 

FUNCTI0N ETA (X) receives the independent variable X = N^ and returns 
the fin effectiveness ETA. The two functions represented by Eqs. 32 a and b 
apply strictly to zero sink temperature, but the effectiveness is insensitive 
to the sink temperature and the equations hold sufficiently well for 0 g < 0.8 0 b * 

ETA (X) is called from MAIN and SDRV. 

f. The FUNCTI0N DETA Subprogram computes the derivative dTj/dN^ of the 
fin effectiveness TjCN^) for the unobstructed fin. The derivative is derived 
from Eqs. 32 a and b, and evaluated by calling the FUNCTI0N P0LY subprogram 

[13. 

FUNCTI0N DETA (X) receives the independent variable X = N c and 
returns DETA = dT/dN^. DETA (X) is called from MAIN and SDRV. 

g. The FUNCTI0N P0LY Subprogram [1] serves to evaluate power poly- 
nomials f(x) of any degree (N - 1). 

N 

f(x) = £ A. X 1 
i = 1 

by computing the recurrence expression 
f i = A N 

f k+l f k X + A N-k ,k = 

(N-l) times. The argument list of FUNCTI0N P0LY (N, A, X) transfers to the 
subprogram the number N of coefficients A_^, the one-dimensional array A(N) and 
the independent variable X; it returns P0LY = f(x). P0LY is called from ETA 
and DETA. 
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3. Input Data Preparation 

The non-dimensional groups 


THETAS =0 = T /T 

(19) 

s so 

U - ’V H Gz 

(22) 

V = (i/n) (t/H) (k/k f ) (N c /N N u ) 

(21) 

FNCBR = N = 2eaT 3 H 2 /(kt) 
c o 

(20) 

RH0D = X » 4H/d 

DXWRT = Az/L, output intervals 

(13) 


are computed and punched, in the above sequence, on two cards in 

F0RMAT (5F16.8) 

There may be arbitrarily many pairs of cards as specified above, 
one pair for each desired simulation. They will be executed in succession. 
The End-of-Job card will cause normal exit from the program. 

The expected run time on the UNIVAC 1108 is 0.24 seconds per 
simulation plus necessary compiling and collecting times. There were 32 
runs executed in 14.8 seconds. 

4. Output Presentation 

For each simulation are printed, on a separate page, first the 
governing parameters as read in: 

THETAS = 0 g , defined by Eqs. 4 & 19 


U 

9 

It 

Eq. 

22 

V 

9 

II 

Eq . 

21 

FNCBAR = 

\ ’ 

II 

Eq. 

20 

LAMBDA = 

X , 

II 

Eq. 

13 


and then a table consisting of seven columns. In the first two columns are 
listed, as functions of axial distance ZETA = ^ at the selected intervals 
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A? = DXWRT, the coolant fluid temperature THETAF = 0^ and the fin base 
temperature THETAB = 0^. The positions ZETA = £ are listed in the sixth 
column. 

The conduction parameter NC = N c , the combined fin-plus-tube 
effectiveness ETABAR = n, defined by Eq. 11, and the fin effectiveness 
(unobstructed fin, Eqs. 32 a and b) ETA = rt, are listed, respectively, in 
the third, fourth and fifth columns. The last columns indicate the inte- 
gration steps required to reach the respective axial position £ . The desired 

coolant fluid exit temperature 9 (U,V,N , X; 0 ) is found in the last entry of 

0 c s 

the first column. 

Figure 11 shows a typical output listing 
C. Minimum Area Optimization 


1. Objective 

The Minimum Area Optimization (MAO) code serves to find the minimum 
coolant fluid exit temperature for a given projected fin panel area, on the 
basis of the optimization analysis developed in Chapter III and the constraints 
as discussed in Section III-2. From the results, evaluated for a number 
of selected areas, one can determine the minimum area required for a 
chosen heat rejection rate. The optimum parameter set X is sought which 
satisfies Eqs. 33, 36, 40 and 41 by performing the iterations specified in 
Eqs. 48 through 53 until Eq. 33 is satisfied. 

From the optimum parameter set one computes the normalized system 
dimensions L*, H*, t* and d* in accordance to Eqs. 59 through 62 or Eqs. 64 
through 67, depending on the flow regime in the coolant channels. 

The necessary derivatives occurring in Eqs. 50 and 51 are obtained 
by Runge-Kutta Integration of the first-order, ordinary, non-linear differ- 
ential equations given by Eqs. 52 and 53. 

The individual program units are discussed in the following sections. 
Reference is made to program description in Section 3 of this Chapter and 
to Reference [1]. 
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2. Deck Assembly 

In its current state the Minimum Area Optimization (MAO) code 

consists of 


1 


main program 

MAIN 


14 

SUBROUTINE 

subprograms 

RKS | 

SDRV . 

SCNTL .! 

I” integration 



* 

BNDCND | 
RESBND | 
REST I 

| restraints 




NDERV 
MDER | 

NCDERV | 

partial 




MIXDER ' 

DERV 

i derivatives 




INITIAL } 

initial 




INTMIX f 

conditions 




FMINV 

matrix inversion 

5 

FUNCTION 

subprograms 

ETA 

DETA 

1 n and its 




DDETA 

DDDETA 

P0LY 

> derivatives 


The basic concept of the program is as described in Chapter B 
before and the block diagram in Figure 7 applies in principle. The 
exceptions are 

(i) there are 20 simultaneous differential equations to be 
solved: 

2 for 0 f and 0 b 

6 for the derivatives of and © b with respect to 
U, V and N c 

12 for the second derivatives of 0^ and 0-, with respect 
to U, V, and N c , D 

(ii) the initial conditions on the last 18 differential equations 
are obtained in subroutines, 
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(iii) partial derivatives occurring in the differential equations 
for the last 18 of the above derivatives are computed by 
subroutines , 

(iv) at the end of each integration, the system of Eq. 49 is 

solved by FMINV and the integration repeated until the optimum 
is reached. 

a. The MAIN Program perforins the same functions as described in 
Chapter IV, Section B.2.a. and, in addition, calls INITIAL and INTMIX which 
compute the remaining 18 initial conditions for the partial derivatives as 
defined through Eq. 54. Moreover, instead of returning to the next input 
data set after completion of an integration, new system parameters X are 
computed by solving, via FMINV, the system of Eq. 49. The results are tested, 
through SUBR0UTINE BNDCND, to satisfy the Inequalities 36. If necessary, 

the newly computed parameters are adjusted to fall within the above limits, 
and the integration is repeated. 

b . The SUBR0UTINES RKS , SDRV and SCNTL perform the integration as 
explained in Sections B.2.b, c and d. There are, however, 18 additional 
dependent variables Y(3) , ..., Y(20) to be integrated by RKS and, consequently, 
18 additional derivatives to be evaluated in accordance with Eqs. 52 and 53. 
These derivatives are evaluated by calling DERV for 30^/3X^, 30^/ 3X^, 3 (i) 2 0^/3X^ 2 
and 8 2 0 b /3X i 2 and by calling MIXDER for 9 2 0 f /9X.,9X^ and B^/ax^Xj, i 1 j. 
SCNTL contains a slightly altered output specification. 

c. Th e Restraint SUBROUTINES BNDCND, RES BNP and REST serve to, 
respectively, 

(i) test newly computed components for compliance with Inequali- 
ties 36. 

(ii) solve Eqs. 40 and 41 for x_ and x. . 

3 4 

When newly computed components X do not comply with Inequalities 36 then 
they are set to meet the minimum and maximum conditions and the system of 
Eq. 49 is solved again, all from within BNDCND. 
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SUBR0UTINE REST is the only particular program which differentiates 
the two optimization modes, minimum area and minimum mass. It may be replaced 
to accommodate other constraints. BNDCND and REST are. called from MAIN, 


d. The SUBR0UTINES NDERV, MDER, NCDER, MIXDER and DERV are used to 
compute , respectively , 

(i) N and its derivatives dN/dA, d 2 N/dA 2 from Eq. 17, 

(ii) M and its derivatives, up to third order, with respect to 

its arguments A and N^ , in accordance to Eqs. 12, 15 and 16, 

(iii) N and its derivatives 3N /30 , 3 2 N /30, 2 and 3 2 N /3T3N , 
c c b c b c c 

(iv) the mixed derivatives defined by Eq. 53 for i ^ j, 

(v) the derivatives defined by Eqs. 52 and 53. 


e. The SUBR0UTINES INITIAL and INTMIX serve to compute the initial 
conditions for the dependent variable Y(3), . .., Y(20), at the channel inlet 
Z = 0 which represent simple and mixed derivatives 3 2 0^/3X^3X^ and 3 2 0^/3}L3X. 
in Eq. 53. 


f. The SUBR0UTINE FMINV is capable of performing two related tasks: 

(i) to invert a square, invertible matrix A 
(ii) to solve a non- trivial system of equations AZ = X. 

The solutions are obtained through a sequence of elementary row 
operations which lead from the properly augmented coefficient matrix to the 
row-reduced echelon matrix, a standard procedure described in most elementary 
introductions to linear algebra [1]. 

SUBROUTINE FMINV (A, X, N, M) accepts, through its argument list, 
the square, invertible matrix A of rank N, and if task (i) above is intended, 
M is set equal to M = 2N. Then, upon return from FMINV there will be the 
inverted matrix A -1 placed in XMAT(I,J) with 1=1, ...,N, J = N + 1, 

N + 2, ..., 2N. The two-dimensional array XMAT is transferred to the calling 
program via a C0MM0N declaration. When task (ii) above is intended, M is 
set equal to M = N + 1 and, upon return from FMINV, the unknown vector Z 
is placed in the one-dimensional array X. 

FMINV is used to solve Eq. 49 and is called from MAIN and from 


BNDCND. 
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g. The FUNCTI0N subprograms ETA» DETA and P0LY are described in 
Sections B.2.e, f, and g of this chapter. 

h. The FUNCTI0N subprograms DDETA and DDDETA are used to compute 

the second and third-order derivatives d 2 n/dN 2 and d 3 n/dN 3 as derived from 

c c 

Eqs. 32 a and b. The derivatives are evaluated via FUNCTION P0LY subprogram. 
3. Input Data Preparation 

Each optimization run is carried out for a selected ASTAR = A 
(Eq. 38) and TSTAR = t* (Eq. 41) as the independent variables. The necessity 
of specifying t* is not yet completely established at this time but included 
in the program for additional flexibility and to prevent the search from 
reading excessively large values of t*. 

The optimization produces 

(i) the optimum coolant exit temperature 
(ii) the optimum parameters X = {x^, ..., x^} from which 
to compute L*, H*, d* and X (Eqs. 59-67). 

For each optimization one computes 


THETAS = 0 

s 

, defined 

by 

Eqs 

. 19 

U, starting value 

II 

3 

if 

Eq. 

22 

V, starting value 

If 

3 

II 

Eq. 

21 

FNCBR = N c 
starting value 

ft 

3 

ft 

Eq . 

20 

ASTAR * A* 

ff 

3 

ft 

Eq. 

40 


and selects 


DXWRT = A? the axial interval for which results are 
to be printed. 

The first five values are punched on one card in 


F0RMAT (5F16.8) . 
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The last value is punched on a second card in 

‘FORMAT (FI 6. 8) 

There may be arbitrarily many pairs of cards as specified above, 
one pair for each optimization run. They will be carried out in succession. 
The End-of-Job card will cause normal exit from the program. 

The expected run time per optimization is approximately 15 seconds 
plus necessary compilation and collection times. 

4. Presentation of Results 

Each iteration produces one page of output, first the three 
parameters THETAS, ASTAR and TSTARI as read in. Next the current values of 
the system parameters U, V, FNCBAR = N^, LAMBDA = X and the identification 
of the flow regime, followed by the iteration' count for establishing initial 
conditions 

The. table lists, as functions of £ - Z/L, in this order 

THETAF = 8f(c) coolant temperature 

THETAB = 0^ (?) fin base temperature 

NC = N (C) = N .6, 3 , Eq. 14 

C CD 

ETABAR = n, defined by Eq. 11 

ETA = P* defined by Eq. 32 

M defined by Eq. 12 

N defined by Eq. 17 

X = £ , the non-dimensional distance along the channel. 

Following the table is a list of the three system parameters 
U, V, and the first and second derivatives of 0^ with respect to these 
parameters. Finally are listed the newly computed changes of the above 
parameters. 

A typical listing of the results is shown in Figure 9. 
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D. Minimum Mass Optimization 

1. Objective 

The objective of the Minimum Mass Optimization . (MMO) code is 
identical to that discussed in Section C above, except that the constraints 
are replaced by the constraints developed in Section III-3 of this report; 
the minimum mass is thought. Equation 46 replaces Eqs. 40 and 41. 

2. Deck Assembly 

The MMO code is assembled identically to the MAO code described 

in Section C.2 above except the SUBR0UTINE REST is replaced by the SUBROUTINE 

C0NST which serves to compute x^ = X from Eq. 46, together with all the 

necessary derivatives of X with respect to x, = U, x„= V and x„ = N . C0NST 

1 2 3 c 

is called from MAIN. RESBND is not required. 


3. Input Data Preparation 

Each optimization is carried out for a selected M ' as the 
independent variable and for the two parameters N* u and > defined by 
Eqs. 46, 42 and 57, respectively. The optimization produces 

(u) the optimum (minimum) coolant exit temperature 
(ii) the optimum system parameters X from which to compute 
(L*, H*, t*, d*} through Eqs. 59-67. 


For each optimization are computed 


THETAS = 0 defined 

s 

^ | starting values 
FNCBAR = N 


c 

T0TMAS = M* , selected 
tot ' 


by Eqs. 


PHI2 = $ 2 
ALPHA = a = 


N* 


Nu 



4 & 19 

22 

21 

20 

46 

44 
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and selected 

DXWRT = At the axial interval for which results are to 
be printed* 

DA, DR absolute and relative integration step errors* 

The data values are punched on the data card in NAMELIST format* 

The NAMELIST is called INPUT. 

There may be arbitrarily many additional input data cards as speci- 
fied above, or with only one datum, one NAMELIST for each optimization run. 
They will be carried out in succession until an End-of-Job card causes normal 
exit from the program. 

The expected run time per optimization is approximately 20 seconds 
plus necessary machine preparation times. 

4* Representation of Results 

Results are presented in the format identical to that described 
in section C.4 above, except that instead of the single TSTAR1 are listed 
MSTAR = M*, PHI2 = <f> 2 and FNUSTR = N* Nu> as read in. The representative 
output example is shown in Figure 10. 
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Figure 13, Results Listing of Minimum Mass Optimization 
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V. CONCLUSIONS 


The work presented herein resulted in a simplified radiator system 
analysis and a systematic optimization procedure. 

Comparison of the simplified with the rigorous [1] system analyses 
indicated that the agreement between the two analyses can be expected to 
be within approximately 5%. 

The optimization procedures carried out lead frequently to an optimum 
on the boundaries of the parameter domain. Mass optimization tends toward 
widley spaced tubes between thin fins. 

Optimization was originally intended to be achieved through the 
parametric study of trends rather than through analytical procedures. 
Within the resources provided by this contract two analytical iterative 
optimization procedures were developed beyond the original goal of work. 
These procedures lead to the maximum heat rejection for given system area 
or system mass. Future work should be aimed toward the development of 

(1) a parameter domain within which relative 
extrema exist, 

(2) suitable working diagrams, through repeated 
applications of the developed codes, which 
depict the optimum geometric system parameters 
and produce either the least area or the least 
weight requirements for a given heat rejection rate. 
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APPENDIX A 

COMPUTER CODE FOR SIMPLIFIED 
RADIATOR SYSTEM SIMULATION 
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APPENDIX B 
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